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Buffon Noodles 


Ed Waymire 


1. INTRODUCTION. During the 1992 spring term, physicist Corinne Manogue 
gave a colloquium talk to our mathematics department in which she surveyed 
various mathematical notions which arise in a string theorist’s approach to quan- 
tum field theory. All this talk of strings and probability by Corinne made an 
impression on my colleague Robby Robson, who had an interesting hallway 
question the following day. Robby was wondering if I knew what would happen if 
Buffon had tossed a noodle in place of a needle (actually Buffon tossed baguettes)? 
Here the image of noodles is that of the long thin stringy type which always seem 
an entangled mess. So this story will feature entangled noodles. 


2. BUFFON NEEDLE. Let’s recall the Buffon needle problem. Snell ([7]) has a 
very nice treatment together with original references and a wonderful dose of 
historical perspective blended with modern computer simulations; that’s where I 
learned about the baguette. One tosses a needle of length L onto a plane surface 
(floor) marked by parallel lines of width D > L units apart. One asks for the 
probability of the event C that the needle intersects its closest line. The answer is 
a (linear) function 


Po (X) ■= P(C) = ^ = y 0 x (1) 

of the ratio x = L/D with slope y 0 = 2/n ( ~ .636). A particularly interesting 
aspect of this example is that one obtains a statistical method of estimating 
numerical values of ir from the crossing frequency of repeated needle tossings. 

One needs a model of the needle tossing experiment to compute the probability 
in (1). For this, let Y denote the vertical distance from the midpoint of the needle 
to the nearest parallel and let 0 denote the acute angle made by the needle and 
the nearest parallel. Then one assumes that (T, 0) is uniformly distributed over 
the range of values [0, D/2] X [0, 7t/2]. Now the problem is a geometry problem. 
It is interesting for the ensuing discussion that the crossing probability (1) does not 
change if one assumes that an endpoint of the needle is randomly (uniformly) 
selected on a fixed vertical line and, independently, the angle made by the needle 
with the vertical line is obtained by a random rotation from 0 to 277 . 

3. BROWNIAN NOODLE. While others (eg., see [1], [5]) have considered “noodle” 
extensions defined by bending needles into various specific convex shapes to be 
tossed, we shall allow the noodles to become randomly tangled . To begin, let us 
first consider what happens in the case of a two-dimensional standard Brownian 
motion {B(t) = (B l (t\B 2 (t))} over a time interval of unit length and started at 
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random along a vertical line of length D , Figure la. Although the total length of 
the Brownian noodle is oo, one may introduce L as a scale parameter by taking 

/ 2L 

5(0 s =— B(t), 0 < t < 1 . ( 2 ) 

The answer in this case is 

pm : = nn = - 2 i - 1 2 (i - i)2(7r2 /4x^/D) 2 ) (3) 

77 m- 1 (2m - 1) 

Now the parameter L is a scale parameter and the gap is D. The crossing 
probability in (3) depends only on the ratio L/D. The extra factor yfl /2 is being 
introduced here for the convenience of making comparisons with computations 
later. For convenience now let L := yflL/2. 



Figure la. Sample Realization of Brownian Noodle Toss. 



Figure lb. Sample Realization of 100 Strand Noodle Toss. 
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To obtain the solution (3) one first observes that it is enough to consider the 
one-dimensional Brownian motion {BjO)} on the interval [0, D/L'] with uniform 
initial distribution on [0, D/L!]. We will need the joint distribution of (M x ,m x ) 
where 


M x = max B x (t ), m x = min B x (t). (4) 

0<f<l 0<f<l 

For then the answer is simply 1 — P ( 0 < m x < M x < D/L). There are various 
ways to compute the joint distribution of the maximum and minimum in (4). One 
approach is to first compute the transition probabilities of the Brownian motion on 
[0, D/L] with absorbing boundaries at 0, D/L and starting at 0 < x <D/L. 
Since these transition probabilities satisfy the heat equation on the interior of the 
interval with Dirichlet boundary conditions one may compute the following eigen- 
function expansion given in ([2], p. 412—413): 


L * ( m 2 ir 2 t I L \ 2 \ ImrrLxX (rmrLy\ 

pO-,x,y)-2-Z^pl — 5 — j 


D 

° < X,y K U I' ^ 


Now the desired probability (3) is the complementary probability to that obtained 
by performing the indicated integrations in 

p|o <m l <M l < — j = —j^ /L f° /L p(l;x,y)dydx. (6) 

Here one also uses the familiar fact that the sum of reciprocal squares is 7t 2 /6, so 
that restricting the sum to reciprocals of squares of odd numbers yields 

7T 2 1 7T 2 7 T 2 

~6 ~ 4 ~6~ “ T* 

Figure 2a contains a plot of the two probabilities (1) and (3) as a function of the 
ratio x = L/D. The corresponding result of a computer simulation is given in 
Figure 2b. 



Figure 2a. Theoretical Brownian Noodle Crossing Probabilities p x (x). 
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Figure 2b. 5000 Toss Simulated Crossing Frequencies for 100 Strand Brownian Noodle. 


4. “REAL” NOODLES. Next let us consider a real string of n needles of lengths 
L/ {n each, and strung together in independent random directions. The total 
length is then n (L/ yfn) = Lyfn . In order to preserve the ratio as x = L/D , the 
gap between parallel lines is chosen as D^n . Since we already know the answer 
for n = 1 by (1), let’s do the problem for n » 1. In this case we first note that the 
functional central limit theorem provides an approximation to the distribution of 
the string of needles by the Brownian noodle. In particular, observe that the 
lengths of the projections of the needle are i.i.d. with, u p to the scale factor 
L/ {n , the arcsine distribution P(Y e dy ) = (1/tt\/ 1 - y 2 )dy. This distribution 
has mean 0 and variance 1/2. Thus the Brownian motion approximation has 
diffusion coefficient L 2 / 2. This explains the scaling factor (v^2 /2)L introduced 
earlier. 

5. FIRST ROUND APPROXIMATIONS. In the limit the approximation involves a 
noodle {SO): 0 < t < 1} of length oo in a gap of width oo, but in the ratio x = L/D. 
This is how string theorists discussing quantum field theory can sound! One might 
“expect” the noodle crossing frequency to be estimated by p x (x/ 4n) as this 
corresponds to formally replacing D by 4nD in (3). This curve is given in Figure 
3a with n = 100. Figure 3b provides the crossing frequencies p 5 (x) obtained from 
a simulation of 5000 tosses with n = 100 for comparison. However, since the effect 
of fixing the gap to length ratio in the above computation is a second rescaling of 
the string of noodles by 1 / 4n , the central limit theorem approximation by the 
Brownian noodle is not justified. On the other hand, the formula (3) is applicable 
to* strings composed of a large number n of needles of (fixed) lengths x each 
tossed onto a surface with gap width 4n . This fact is illustrated by the simulation 
in Figure 2b of the corresponding crossing frequencies pfx) for a Brownian 
noodle with n = 100, N = 5000; i.e. these are the crossing frequencies of - 4n /2 
or yfn / 2 by a string of n randomly (uniformly) oriented needles of length x each 
with endpoint randomly distributed over [ - 4n /2, 4n /!]. 

Due to the additional rescaling of the gap width, there will be no loss in 
generality to consider the crossing probability of strings composed of n needles of 
length x < 1 each and a gap of width n. A quick and clean (rigorous) upper bound 


1994] 


BUFFON NOODLES 


553 




Figure 3a. Theoretical Brownian Noodle Crossing Probabilities p x (x/ yf\ 00). 



on the crossing probability can be obtained with Doob's maximal inequality ([2], 
pp. 49-52). In particular one obtains, using translation invariance of the problem 
to pose it with one end started uniformly in [ — n/ 2, n/ 2], that 

4 o 2 

P(C) < -x 2 + -. (7) 

For large n this bounds the probability by the essentially constant value 2/3+ . 
Recall that the maximum crossing probability for the needle (1) is 2/ir ~ 2/3 - . 
To make this estimate, let {S n } denote the random walk starting at 0 with 
displacement distribution P(Y e dy) = ( X/iryjx - y 2 )dy , and let Z (n) = nU be 
uniformly distributed on [-n/2, n/2] independently of the random walk. The 
details of the distribution of Y do not play a role beyond first and second 
moments. The essential observation is that by symmetry {xS k + Z (n) : k = 0, 1, ... } 
is a martingale. Thus, Doob’s inequality gives (7), after conditioning on Z (n) and 
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using symmetry, as follows: 


P(C) = £p|max({xS* + Z (n) ) >^j U jminfxS*. + Z (n) ) < -^j|Z (n) j 

< EP [max (xS k + Z w ) > ^|Z (n) ) + Ep[ min (xS k + Z (n) ) < - ^| Z (n) ) 
= 2£/ > (max(xS yt + Z< n) ) ^ ^|Z (n) ) < 2^£(xS„ + Z (n) ) 2 

( 8 ) 


_ 

8 1 

x 2 n 

' 

h £) 

-i x 2 + 

2 

! 



12 ] 

n 
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A bound as universal as (7) will be much too large for many cases of interest. More 
precise estimates of the crossing probability upper bound are given in the next 
section. 

For a lower bound, let ZX/3) denote the event that the initial end of the string 
falls within /3x, /3 < 1, units from the boundary. Then D(fi) has probability 
2fix/n. Also, from the geometry P(C |D(/3)) > (Cos _1 (/?)/7 r) since, for any n, 
crossing is implied by the occurence of the angle of the first needle between 
+ Cos -1 (/3). Thus, for any 0 < /3 < 1, P(C ) > 2/3x/rt(Cos _1 (/3)/7r) so that 

2 i .357 

P(C) > \{n)x, where A (n) = — max /3Cos J (/3) ~ . (9) 

/1 77 0 </3 < 1 /I 


6. LARGE DEVIATIONS. As noted above, the functional central limit theorem 
approximation involved in replacing the string of needles by a Brownian motion for 
large n provides probabilities of events involving 0(1) fluctuations in the shape of 
the string of needles; the individual needle lengths already being 0(1/ v^). 
However, the crossing probabilities of 4nD involve 0{yfn) deviations from this 
picture which are not obtained by this approximation. These fluctuations have 
probabilities which also depend on the details of the angular distribution of the 
needles in the string. The precise probability computation desired here is a large 
deviation problem of the type which arises in the actuarial mathematics of insur- 
ance risk. In fact, the classic work of Cramer on large deviations was motivated by 
actuarial problems. Posed this way, the noodle problem also has an interesting 
twist in that it involves a “two-sided ruin event.” 

To see how the distribution will depend on the angle distribution, consider the 
degenerate case in which the angles made by the needles comprising the string 
with the horizontal are each ir/2 (nonrandom). Then the noodle is a needle with 
fixed orientation and the only randomness is in the distribution of the endpoint; 
see Figure 4a. For this case one readily finds that p 3 (x) = P(C ) = x . (This 
prgbability is unchanged if one randomizes the choice of endpoint). To see how 
large deviation estimates go in a simple but less trivial situation for which explicit 
computations are still possible, consider the case in which the needles comprising 
the string make angles ±7t/4 with the horizontal with equal probabilities; see 
Figure 4b for a sample outcome of a simulation of a single toss with n = 100. 
Since the variance in step size is one in this case, the rescaling by ^2 /2 is 
unnecessary for the comparison with the (nonrigorous) central limit theorem 
prediction; i.e. replace x by x/2 in Figure 3a to correct for the ^2/2 before 
comparing with the observed Figure 5. In this case the lengths of the projections 
of the needle are i.i.d. with, up to the scale factor L, a Bernoulli distribution 
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Figure 4a. Sample 0 Degree Single Strand Noodle Realization. 



Figure 4b. Sample 100 Strand 45 Degree Noodle Realization. 


P(Y, = + 1) = P(Y t = - 1) = 1/2. Let S n ■■= (Y, + Y 2 + ■■■ + Y„), S 0 = 0. More 
generally, assume (f>(s) = Ee sY ' finite and note for i.i.d. { Y„}, (<£(s))” = Ee sS " > 
e" sb P(S n > nb). That is, for arbitrary s 

P(S„ > nb) < € "('oh(^))-^. (10) 

In particular, 

P(S n > nb) < e" x * (b) , where x*(b) = inf(log Ee sV - sb). (11) 

S 

The reverse (asymptotic) inequality for b > EY X yielding 
1 

— log P(S n > nb) -> x*(b) as n -+ <*>, (12) 


is the often cited Cramer-Chemoff theorem on large deviations of S n /n from EY V 
That is to say, P(S n > nb) ~ e nx * {b) as n -> oo; the rate ~x*(b) is a large deviation 
rate. The reader may consult the developments given in ([3], pp. 147-148) and ([4], 
pp. 57-63) for some basics of large deviation theory. 

Using the reflection principle for simple symmetric random walk one may also 
check that ([2], p. 10, 70), PGnax^,, S k > b) = 2 P(S n > b) - P(S n = b). In par- 
ticular one obtains from this, 

p( maxS k >b) < 2 P(S n > b). (13) 

V k<n } 
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These results will serve as our basic tools. Notice that for large n it is no more than 
a venal sin to think of the inequality (13) as equality; i.e. the resulting probability 
bound may be expected to be close to the actual probability. 

First let us compute the so-called Legendre transform x*(b) of the function 
*(s) = log Ee sY required in (11). One may easily check by calculating limits as 
s -> ±oo that x*(b) = -oo if \b\ > 1. For \b\ <1, equating the derivative with 
respect to 5 of {log((e 5 + e~ s )/2) - sb } to 0 and solving for 5 gives e 2s = ((1 + 
b)/( 1 - b)) and therefore 

X*( b ) = ~ -y-log(l + b) - !=^log(l - b ). (14) 

As before let Z (n) = nU be uniformly distributed over [~n/ 2, n/2] and indepen- 
dent of {S n }. Then one has along the lines of (8) that 

/.(C) £ + 7 }) s 4f(s» i »{i + 7 }). (15) 

first using the symmetry of the distribution of both {.S',,} and U, and then the 
reflection principle. By conditioning on U, using (15), the argument for (11) 
conditionally, and then (14), one has after taking expected values (with respect to 
U), 

P(C) < AEP ( S„ > n(L + ^J|i/J < 4 Ee nx*u/ 2 x+u/x) _ 4 X j 1/X e nx * w dy 


4*J^expj--[(l + y)log(l + y) + (1 -y)log(l -y)]|dy : = y 4 (n)x. 

(16) 


Note that the (concave) function 

*(y) = — ( 1 / 2) [(1 +y)log(l +y) + (1 -y)log(l - y)] 

and its first derivative are zero at y = 0. The second derivative at 0 is —1. As a 
result one may bound by a Taylor’s approximation -(l/2)y 2 to obtain by a simple 
change of variable and the fact that 

f^ e -o/2)P dz „ ^ 2 ^/ 2 , 

• / o 

y 4 (n) = 4^ 1 exp|- -[(1 + y)log(l + y) + (1 - y)log(l - y)]| dy 

< 2\l2^n~ l/1 * 5.01«- 1/2 . (17) 

Observe from this that the crossing probabilities are theoretically predicted to be 
smaller than the universal bound (7) for 45 Degree strings consisting of n strands. 
FiduRE 5 depicts the crossing frequencies p 4 (x) for a simulation of N = 1000 
tosses with n = 100. A sample realization of a toss was plotted in Figure 4b. 
Finally, observe that the argument used to arrive at the lower bound (9) can be 
simply adapted to this case to give 

x 2 v^tt 

- <P(C) < y 4 (n)x -j^x. (18) 

n V« 

To extend these estimates to the case of uniform angles on [0,2v) first requires 
a general inequality of the form (15). Such an extension is possible for symmetric 
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Figure 5. 5000 Toss Simulation Crossing Frequencies for 100 Strand 45 Deg Noodle. 


random variables ( Levy's inequality ), ([5], p. 14; [2], Exc. #4, p. 70). The corre- 
sponding bound in place of (16) for this case is, following the first two inequalities 
of (16), 

P(C) < 4 Ee nx * (l/2x + u/x) = 4x [ 1/X e n ** (y) dy = 4x[ 1 e n ** M dy, (19) 

J o J o 

where now for |y| > l,**(y) = — oo, but for |y| < 1 



The crossing frequencies p 5 (x) obtained from a simulation of 5000 tosses with 
n = 100 strands were provided earlier in Figure 3b. To obtain the asymptotic 
value of y 5 (n) one may proceed as in the above case to compute the values of 
**(y) and its first two derivatives at y = 0. This is most easily accomplished by the 
following Legendre transform duality equations : Let 



Then 

dx{s) dx*{y) 

**0) = x(s(y)) -ys(y), — = y, — — = -s. (22) 

One may check from this that **(y) and its first derivative are 0 at y = 0 and the 
second derivative is - 1 /2. It follows as above that 

y 5 (n) =4 f e nx * (y) dy < 4yfrrn~ 1/2 ~ 7.09n _1/2 . (23) 

•'o 

A lower bound was given at (9). In this case the crossing probability prediction falls 
below the bound (7). 

What we have presented here is the stuff with which probability bounds and 
estimates can be made. The computer simulations provide an order of magnitude 
feel for crossings, but no more than this for n large. The reader desiring a more 
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precise recipe for “n-noodle pi” may want to give more careful consideration to 
the small n cases, beginning with n = 2. We only hope to have whet the appetite 
for more probability and noodles. 
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PICTURE PUZZLE 

(from the collection of Paul Halmos) 



What can these men talk about? 
(see page 570) 
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